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Magnetization and its relaxation have been measured in Ba(Fei_i,Coa:)2As2 single crystals at var- 
ious doping levels ranging from very underdoped to very overdoped regimes. Sizable magnetization 
relaxation rate has been observed in all samples, indicating a moderate vortex motion and rela- 
tively small characteristic pinning energy. Detailed analysis leads to the following conclusions: (1) 
A prominent second-peak (SP) effect was observed in the samples around the optimal doping level 
{x ~ 0.08), but it becomes invisible or very weak in the very underdoped and overdoped samples; 
(2) The magnetization relaxation rate is inversely related to the transient superconducting current 
density revealing the non-monotonic field and temperature dependence through the SP region; (3) 
A very sharp magnetization peak was observed near zero field which corresponds to a much reduced 
relaxation rate; (4) A weak temperature dependence of relaxation rate or a plateau was found in the 
intermediate temperature region. Together with the treatment of the Generalized-Inversion-Scheme, 
we suggest that the vortex dynamics is describable by the collective pinning model. Finally, vortex 
phase diagrams were drawn for all the samples showing a systematic evolution of vortex dynamics. 

PACS numbers: 74.25.Qt, 74.25.Ha, 74.70.Dd 



I. INTRODUCTION 

Since the discovery of superconductivity at Tc — 26 
in LaFeAsOi_a:Fa;, the iron-based layered superconduc- 
tors have exposed an interesting research area to the com- 
munity of condensed matter physics. This breakthrough 
was followed by the realization of superconductivity at 
temperatures at as high as 56-57 K in Cai_2:R.EFeAsF 
(fluorine derivative 1111 system with RE = Nd, Sm, 
etc.).— A lot of experimental and theoretical works on the 
physical properties have been carried out. For the FeAs- 
1111 phase, it is very difficult to grow crystals with large 
sizes, therefore most of the measurements so far were 
made on polycrystalline samples. This is much improved 
in the (Ba,Sr)i_2;Ka:Fe2As2 and (Ba,Sr)(Fei_2;Coa;)2As2 
(denoted as FeAs-122) system since sizable single crystals 
can be producedi^ Specific heat, lower critical field, mi- 
crowave measurements suggest that these superconduc- 
tors exhibit multiband featurei'^'S. Measurements un- 
der high magnetic fields reveal that the iron-based su- 
perconductors have very high upper critical fieldsj^'i'^ 
which indicate encouraging potential applications. Pre- 
liminary experimental results also indicate that the vor- 
tex dynamics in iron pnictide may be understood with 
the model of thermally activated flux motion within the 
scheme of collective vortex pinmngi^ i^°i^^i^^ A second- 
peak (SP) effect has been observed in Bai_3;K2,Fe2As2iS 
and Ba(Fei_a;Coa;)2As2 single crystals Jii^ii^ii^ More- 
over this system has a layered structure with the con- 
ducting FeAs layers being responsible for the supercon- 
ductivity, which behaves in a similar way as the case of 
cuprates. For cuprate superconductors, due to the high 
anisotropy, short coherence length, and high operation 
temperature, the vortex motion and fluctuation are quite 
strong. ""^^ This leads to a small characteristic pinning en- 



ergy, and the single vortex or vortex bundles are pinned 
collectively by many small pinning centers. Therefore 
it is curious to know whether the vortex properties and 
phase diagram of the cuprate and FeAs-based supercon- 
ductors are similar to each other or notJ^'i^'i^ In this pa- 
per, we report an intensive study on the vortex dynamics 
of Ba(Fei_a;Coa;)2As2 single crystalline samples ranging 
from very underdoped to very overdoped regime. 



II. EXPERIMENT 

The single crystals with high quality measured in this 
paper were prepared by the self-flux methodi^ Sam- 
ples with six different doping concentrations {x — 0.06, 
0.07, 0.08, 0.1, 0.12, 0.15) with typical dimensions of 
1.0 X 1.0 X 0.3 mm^ have been used for both magnetic and 
resistive measurements. The measurements were carried 
out on a physical property measurement system (PPMS, 
Quantum Design) with the magnetic fleld up to 9 T. Dur- 
ing the measurements, the magnetic field H is always 
parallel to c-axis of single crystals. The temperature sta- 
bilization of the measuring system is better than 0.01 
K. The magnetic properties were measured by the sensi- 
tive vibrating sample magnetometer at the vibrating fre- 
quency of 40 Hz with the resolution better than 1 x 10~^ 
emu. The magnetic field sweeping rate can be varied from 
0.5 Oe/s to 627 Oe/s, for most of the measurements of 
dynamical magnetization relaxation measurements (de- 
fined below) we adopted the sweeping rate of 50 Oe/s 
and 200 Oe/s, which were much enough for us to resolve 
the difference of magnetization. The advantage of this 
technique is that the speed of data acquisition is very 
fast with a quite good resolution for magnetization. 

In Fig.l we show the temperature dependence of the 



2 




15 30 45 

T(K) 

FIG. 1: (Color online) Temperature dependence of the super- 
conducting diamagnetic moment (upper panel) and resistiv- 
ity (bottom panel) of six doped samples. The inset of upper 
panel shows the temperature dependence of the diamagnetic 
moment measured in the ZFC and FC processes at a field 
of 20 Oe for the sample x — 0.08. The sharp magnetization 
and resistivity transition curves assure the high quality of the 
samples. 



diamagnetic moment and resistive transitions of the six 
doped samples. The superconducting transition temper- 
ature Tc is found to shift systematically with doping, re- 
vealing a dome-like doping dependence (shown later in 
Fig. 20). The sample with nominal composition x = 0.08 
was found to be optimally doped with the highest onset 
transition temperature « 2A.5K. In the underdoped 
region {x < 0.08), an upturn in the resistivity curve above 
Tc can be easily seen, which was supposed to be related 
to the structural and antiferromagnetic (AF) transition. 
The large difference between zero field cooling (ZFC) and 
field cooling (FC) magnetization (shown in inset of Fig. 1) 
indicates a strong magnetization hysteresis in the sample. 
The perfect diamagnetism in the low temperature region 
and sharp transitions observed from the ZFC curves in- 
dicate the high quality of our samples. 

To investigate the vortex dynamics, in this work we 
adopt the so-called dynamical relaxation method i^^'^^i^^ 
Dynamic magnetization-relaxation measurements were 
carried out in the following way: The sample is cooled 
down to a chosen temperature in ZFC mode and then 
we sweep the magnetic field and measure the magnetic 
moment by following the routes: — * Hmax — > — > 
—Hrnax — > — > Hmax with different field sweeping rates 



dB/dt. The corresponding magnetization-relaxation rate 
Q is defined as 

^_ din J, _ dln(AM) 

^ d\n{dB/dt) dln(dS/dt)' ^ ' 

where js is the transient superconducting current density, 
dB/dt is the field sweep rate. Comparing to the con- 
ventional magnetization-relaxation measurement (fixing 
the magnetic field and measuring the time dependence 
of magnetization) , this dynamical relaxation method can 
overcome the following drawbacks: (1) For the samples 
with a large demagnetization factor, a slight overshoot of 
the field (even lower than 1 mT) can modify the current 
distribution dramatically; (2) A long waiting time is nec- 
essary before meaningful relaxation data points can be 
recorded; (3) To get the valid relaxation data, we need 
to measure the magnetization relaxation in a long time 
period. 

In Fig. 2 we show the magnetization hysteresis loops 
(MHL) of six doped samples measured at 2 K with the 
magnetic sweeping rates of 50 Oe/s and 200 Oe/s, re- 
spectively. The symmetric MHL curves indicate that 
the bulk superconducting current instead of the surface 
shielding current dominates in the samples during the 
measurements. A surprising observation here is that the 
difference between M measured at 200 and 50 Oe/s can 
be easily distinguished, which indicates a relatively large 
vortex creep rate, or as called as giant vortex creep in the 
cuprate superconductors. A prominent SP effect or called 
as the fish-tail effect was observed in the samples around 
the optimal doping level {x = 0.08), but it becomes invis- 
ible in the very underdoped sample(x = 0.06) and hardly 
visible in the highly overdoped samples {x = 0.15). When 
the field is approaching zero, the absolute value of mag- 
netic moment increases markedly. The central peaks are 
surprisingly sharp in all the six doped samples which was 
hardly seen in the conventional superconductors and high 
temperature superconductors (HTSC). Whereas in FeAs 
based superconductors, this kind of sharp peaks are of- 
ten observed in both 1111 and 122 systems, which will 
be detailed in next section.— 



III. RESULTS AND DATA ANALYSIS 
A. The optimally doped sample 

1. The sharp central peak and second-peak 

The MHLs measured at different temperatures from 
2 K to 20 K are presented in Fig. 3 for the optimally 
doped sample {x = 0.08). The symmetric curves sug- 
gest that the bulk pinning instead of the surface barrier 
dominates in the sample. The SP effect can be easily 
observed in Fig. 3(b) at T > 10 K. With the decreasing 
of temperature, the SP moves to higher magnetic field 
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FIG. 2: (Color online) Magnetic hysteresis loops M vs H of the six doped samples at 2 K. All the samples were measured 
with the field sweeping rate of 50 Oe/s (shown by red circles) and 200 Oe/s (shown by black squares). One can see that the 
second-peak or fish-tail effect appears near the optimal doping, while it becomes invisible in the very underdoped and very 
overdoped samples. 



and finally goes beyond the maximum field value 9 T as 
shown in Fig. 3 (a). The global shape of MHL and the 
related features resemble that in the cuprate supercon- 
ductor YBaaCuaOr-a.^ 

When the external field is approaching zero, the abso- 
lute value of magnetic moment increases markedly and 
reaches a maximum, showing a very sharp magnetiza- 
tion peak near zero field. In high temperature supercon- 
ductors (HTSC) and other conventional type-II super- 
conductors, we can see a magnetization hump near zero 
field. This hump near zero field was illustrated by math- 
ematical simulation as due to the following reason:^'* The 
low flow rate or strong pinning in the interior part pre- 
vents the vortices from leaving the sample rapidly, while 
when the external field is approaching or just departing 
from zero, the vortices of opposite signs enter the sample, 
which makes the number of vortices decrease quickly at 
the edge. This produces a non-linear B{x) profile with 
a high slope near the edge. Clearly, in our six doped 
samples the central magnetization peaks are surprisingly 



sharp. This sharp magnetization peak near zero field 
in iron pnictide may be understood in the similar way. 
When the external field is swept back to zero, because of 
the small absolute value of B{x), the slope of B{x) and 
the critical current density near the edge is much larger 
than that in the interior part. In this much en- 

hanced magnetization will appear. However, this simple 
picture will be complicated by the step of B{x) in the 
surface layer with thickness of penetration depth. Mean- 
while the surface barrier and/or the geometrical barrier 
may also play roles here. 

Based on the Bean critical state model)^ we calculate 
the superconducting current density 



where AM = M+ - M_, and M+ (M_) is the mag- 
netization associated with increasing (decreasing) field; 
w, I is the width and length of the sample separately. In 
Fig. 3 and Fig. 4, we present the field dependence of js and 
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FIG. 3: (Color online) The MHLs of the optimally doped 
sample measured at 2 K, 4 K, 8 K, 10 K, 15 K and 20 K. The 
solid line was measured at the magnetic field sweeping rate 
of 200 Oe/s while the dotted line at 50 Oe/s. A second-peak 
effect appears in all the MHLs. 



relaxation rate Q. As shown in Fig. 3 and Fig. 4, the mag- 
netization relaxation rate is inversely related to the js re- 
vealing the non-monotonic field and temperature depen- 
dence through the SP region. Therefore the SP effect is 
induced by the transient relaxation effect which depends 
on the time scalei22i2^ As mentioned above, near zero 
field, a clear sharp magnetization peak can be found in 
MHL. Accompanying this peak, there is a clear suppres- 
sion of the relaxation rate which exhibits a valley near 
zero field. This may be related to the stronger critical 
current density in the edge region where the B{x) value 
is small. With increasing field, a peak of magnetization 
relaxation can also be observed region corresponding to 
the crossover of js between the low-field high slope and 
high-field low slope. This may suggest that the crossover 
point (or the peak position of Q) corresponds to two dif- 
ferent regimes of vortex dynamics. The position of the 
peak shift to the lower field with the increasing the tem- 
perature. A full understanding to this effect would need a 
local measurement facility, such as Hall probe array and 
magneto-optics. We leave this to a future investigation. 

Another interesting feature revealed by Fig. 4 and Fig. 5 
is that, the magnetization relaxation rate rises up quickly 
and monotonically to 100% when the magnetic field is 
beyond the SP position Hsp- This may suggest that the 
high field region is dominated by the plastic motion of 
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FIG. 4: (Color online) (Upper panel) Magnetic field depen- 
dence of the MHL width AM and (Lower panel) dynamic 
magnetization relaxation rate from 2 K to 9 K with the sam- 
ple of a; = 0.08. The SP position corresponds to a minimum of 
the magnetization relaxation rate Q. Near zero field, a sharp 
drop of relaxation rate was observed in accompanying with 
the sharp central peak of MHL. 



vortices. This conclusion will be corroborated by the 
detailed analysis based on the collective pinning model 
in next section. 



2. Analysis based on the vortex collective pinning model 

The temperature dependence of the transient super- 
conducting current density js calculated through AAf 
based on the Bean critical state model and the dynami- 
cal magnetization relaxation rate are presented in Fig. 6. 
ft is clear that the logjs — T curve shows a crossing fea- 
ture at about 10 K for the data measured with different 
magnetic fields, this is understandable since the SP effect 
exhibits a crossover between different pinning regimes. 
Interestingly the relaxation rate Q shows a weak tem- 
perature dependent or a bump-like behavior in the in- 
termediate temperature region. Although the relaxation 
rate Q exhibits a rather large value in the low temper- 
ature approach, we cannot conclude whether the quan- 
tum tunneling of vortices is strong or not in the pnictide 
superconductors since the lowest temperature measured 
here is about 2 K. The plateau or bump- like temperature 
dependence of Q was also observed in cuprate supercon- 
ductors, which is especially pronounced in YBa2Cu307_5 
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FIG. 5: (Color online) Magnetic field dependence of the tran- 
sient superconductivity current js and (Lower panel) dynamic 
magnetization relaxation rate from 10 K to 22 K with the in- 
crements of 2 K for the sample of a: = 0.08. The SP position 
corresponds to a minimum of the magnetization relaxation 
rate Q. Near zero field, a sharp drop of relaxation rate was ob- 
served in accompanying with the sharp central peak of MHL. 
A quick and monotonic rising of the relaxation rate after the 
SP position was observed. 
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FIG. 6; (Color online) Temperature dependence of log jj, and 
dynamic magnetization relaxation rate Q with x = 0.08 at 
0.5, 1, 3, 5 and 7 Tesla. A plateau (or a weak bump) appears 
in the intermediate temperature region. 



superconductori^ This plateau cannot be understood 
within the picture of single vortex creep with the rigid 
hopping length as predicted by the Anderson-Kim model, 
but was attributed to the effect of collective pinning. We 
will illustrate this point in the following discussion. In the 
high temperature region, the relaxation rate rises sharply 
corresponding to the plastic motion of vortices. 

In order to understand the vortex motion in a detailed 
way, we adopt the model of thermally activated flux mo- 
tion (TAFM), which reads 



E — vqB exp(- 



knT 



(3) 
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here E is the electric field induced by the vortex motion, 
U{js, T, Be) is the activation energy, vq is the attempting 
hopping velocity, i?e is the actual local magnetic induc- 
tion. It was suggested that the activation energy can be 
written in a very general way as 



Uijs,T,B, 



^c(r,i?e) .. Jc(T,j?e) .^(T.B. 
/i(r,Be) ^^js{T,B,)' 



1], (4) 



FIG. 7; (Color online) Temperature dependence of the ratio 
T/Q{T) for the optimally doped sample for different fields. 
The slope of T/Q is positive in the low and intermediate 
temperature region, indicative of an elastic vortex motion(see 
next). While it becomes negative in the high temperature 
region revealing the plastic motion of vortices. 
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where fj,, Uc, and jc are the glassy exponent, intrinsic 
characteristic pinning energy and the un-relaxed critical 
current density, respectively. The latter two parameters 
Uc and jc strongly depend on the pinning details, such 
as the characteristics of the pinning centers, the disor- 
der landscape, the condensation energy, the coherence 
length and the anisotropy, etc.. The glassy exponent fj, 
gives influence on the current dependence of U, which is 
a decreasing function of js for all possible values of /z. 
From the elastic manifold theory, it was predicted that fi 
= 1/7, 3/2, 7/9 for the single vortex, small bundles and 
large bundles of vortex motion with the weak collective 
pinning centers. '^'^ For /i = -1 the equation recovers the 
Kim- Anderson model, and the Zeldov logarithmic model 
as special cases can be described for ^=0.— Just by com- 
bining above two equations, it was derived that^ 



^d^ \nE 
' dlnj2 ■ 



(5) 



This equation indicates that a negative ^ would corre- 
spond to a positive curvature of d In i? vs. d In js , indicat- 
ing of a flnite dissipation in the small current limit and 
plastic vortex motion, while a positive /i corresponds to 
a negative curvature of din i? vs dlnj,,, showing a van- 
ishing dissipation in the small current limit and elastic 
vortex motion. Therefore Eq.5 is physically meaningful 
for any value of /i, both positive and negative, as well as 
0. Eq.5 is most suitable for a consistent analysis when 
a transition occurs in the type of pinning as reported in 
this paper. From the general vortex motion equations (3) 
and (4) mentioned above and the definition of Q , Wen 
et al. derived the following cquation^^ 



Q{T,Bc) 



Uc{T,Bc) 



fi{T, Be)CT, 



(6) 



where C = ln{2voB /IdBc/dt) is a parameter that is 
weakly temperature dependent, I is the lateral dimen- 
sion of the sample. We thus present the T/Q vs. T 
at different magnetic fields in Fig. 7. It is clear that the 
curve T/Q vs. T gives a positive curvature in the low and 
intermediate temperature region, while a negative slope 
appears in the high temperature region. Upon above dis- 
cussion, the positive slope of T/Q vs. T would suggest 
an positive glassy exponent fi and elastic vortex motion, 
and a negative slope of T/Q vs. T may suggest a neg- 
ative ^ provided the temperature dependence of Uc{T) 
is not strong. By extrapolating the curve T/Q down to 
zero temperature, one can obtain the value of [/c(0). The 
value of Uc{0)/kB at 0.5 T calculated from Fig. 7 is about 
98 K, which is actually a small value, implying a quite 
small characteristic pinning energy. The Uc{Q) is about 
300 K (at 0.5 T) in YBCO thin mmsM^ but beyond 
3000 K in MgB2.~ With the increase of magnetic field, 
the Uc{0)/kB gradually decreases. While there is an up- 
turn at 7 Tesla, where the SP effect sets in. By calculat- 
ing the slope of T/Q vs. T curve, we can get the value of 
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FIG. 8: (Color online) Temperature dependence of the un- 
relaxed critical current density jc{T,Be) and the normalized 
intrinsic pinning energy g{T, Be) curves obtained by means of 
the GIS, for the optimal doped sample. 



fiC assuming that Uc(T) is not a strongly temperature 
dependent function. In the intermediate temperature re- 
gion, C can be derived from the formula— 



dlnj\i 
dT 



dT 



C— . 

T 



(7) 



Here we found C — 28.56 ±0.66. In this way, fj, can 
also be roughly estimated. At 0.5 Tesla, we obtained a 
positive value /i = 0.45, which indicates an elastic vortex 
motion. The value obtained here by the rough estimate 
is also consistent with what we get through a quantita- 
tive analysis based on the Generalized-Inversion-Scheme 
(GIS). 



3. Generalized-Inversion-Scheme 

To extract information on the unrelaxed critical cur- 
rent density jc{T,Be) and the corresponding character- 
istic pinning energy Uc(T, Be) directly from the relax- 
ation data, the GIS was proposed by Schnack et alj^ 
and Wen et ali^^. The basic assumptions of GIS are: 
(1) TAFM [eq.(l)]; (2)t/(j„r) = UciO)git)f{js/jc); 



7 



(3)5(0 « [jc{T)]PG{T), here UM9{js/jc)G{t) andpde- 
pend on the specific pinning models (see below). These 
assumptions represent the most general scheme among 
many methods proposed so far. According to GIS, one 
can determine jc(T) by using the following integral: 

j,(T) = j,(0) X exp[ 

° CQ(r')[l -dlnG(r')/dlnT'] +dlnjs(T')/dlnr' 



l+pQ{T')C 



dT' 



(8) 



Here jc(0) is the critical current density at K. In this 
equation, the Q(T) and js{T) are the measured values. 
Once jc{T) is known, the temperature dependence of the 
intrinsic pinning energy Uc{T) is obtained from assump- 
tion (3) of GIS. To apply this procedure, one must know 
the value of p as well as the function G{T). For a 2D 
pancake system, it is known that Uc{T) = jciT)(/)Qdrp, 
where (j)o is the flux quanta, d is the thickness of the 
pancake vortex, rp is the pinning range, which typically 
is identified with the coherence le ngth ^. Thus, in this 
special case, p = I and G(T) = ^(1 + t^)/{l - t^) with 
t = T/Tc- Similarly, for a 3D single vortex, it is found^^ 
that p = 0.5 and G(T) = (1 + ^2)5/4/(1 _ ^2)1/4^ gy 
inserting these expressions into above equation, one can 
calculate jc{T) and Uc{T) from the experimental data. 
The results of such a GIS analysis based on the single 
vortex approach is presented in Fig. 8(a) and (b). To be 
consistent with the analysis given in the context of Fig. 7, 
the same value of C— 28 has been assumed. In Fig. 8(b), 
we show together the theoretical predictions for the two 
basic pinning mechanism: pinning due to the spatial fluc- 
tuation of superconducting transition temperatures Tc 
which is called fJTc-pinning,"^^ and that due to the spa- 
tial fluctuation of the mean-free-path, the so-called Sl- 
pinning»i^i2i Taking p = 0.5 and C — 28, the unrelaxed 
critical current density jc{T) and function g{T) are de- 
termined and presented in Fig. 8. Shown together with 
g{T) are the theoretical predictions for the single vor- 
tex of (5Z-pinmng with g{T) ~ 1 ~ and for the 6Tc- 
pinning g{T) — {1 — t'^)'s (1 + t'^)2 , One can see that 
the experimentally derived value resides in between the 
JTc-pinning and (5/-pinning. However an enhancement of 
g{t) in the high temperature region is hardly achieved 
by the JZ-pinning, which is however anticipated by the 
model of JTc-pinning. We would therefore conclude that 
either (5Tc-pinning or some other pinning mechanism are 
in functioning in the present optimally doped samples. 
This of course warrants further clarification. 

Fig. 9 shows the U{js,T, Be) obtained by GIS at five 
fields for the optimally doped sample. The data at 0.5 
T was fitted to eq.4. The fitting parameter Uc is 116 K, 
which is close to 98 K mentioned above. From the fit 
we get /i = 0.55, which is quite close to that obtained 
before fj, = 0.46. We should mention that the fj, value 
determined here is just an averaged one, which in prin- 
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FIG. 9: (Color online) The U{js,T,Be) relation obtained by 
means of the GIS with p = 0.5, for the optimally doped sam- 
ple. The solid line gives a theoretical fit to the U (js ,T,Be = 
0.5r) as formulated by eq.3 with Uc ~ 116 K and fi — 0.55. 
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FIG. 10: (Color online) The vortex phase diagram of the 
optimally doped sample with x = 0.08. 



ciple is also current dependent. Our work suggests that 
the collective pinning model is applicable in this kind of 
superconductors at low and intermediate temperature re- 
gion with a positive glassy exponent. At high fields (still 
below Hgp), the collective pinning model may still work, 
but it is difficult to be quantitatively described by the 
GIS since crossover from the single vortex creep to small 
bundles or large bundles have occured. 



4- The vortex phase 



In Fig. 10, we present the vortex phase diagram of 
the optimally doped sample with x= 0.08. From the 
magnetic measurement, three characteristic fields are de- 
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FIG. 11: (Color online) The MHLs ol the overdoped sample 
{x — 0.12) from 2 K to 16 K. The solid line represents the 
data measured at the magnetic field sweeping rate of 200 Oe/s 
while dotted line at 50 Oe/s. 



termined as shown by the solid symbols in Fig. 10. The 
second-peak Hgp and the relaxation rate peak Hqpeak to- 
gether with the irreversibihty field Hirr (taken with the 
criterion of 0.1 emu/cm^) are shown. The upper criti- 
cal field Hc2 with 95 percent p„ is shown by the filled 
circles. There is a large area between the Hsp — T, and 
Hirr — T curves, suggesting that the vortex dissipation 
is through a plastic motion in this region, but the dissi- 
pation level is still quite low. The Hsp — T, Hgpeak ~ T 
and Hirr — T arc clearly temperature dependent. As an 
accumulation of knowledge, the three curves are well fit- 
ted by the expressions H,p{T) = HspiO)il - T/Tc)i ^ 
HQpeakiT) = HQpeakmi " T/T,)^-^ and H,rr{T) = 
H,rr{0){l - T/Tc)i-5^ respectively. 



B. Overdoped sample with x=0.12 

Now we turn to the case of an overdoped sample {x = 
0.12) with superconducting transition temperature of 19 
K. The MHLs of the sample x = 0.12 measured at 2-16 K 
with the magnetic sweeping rate 50 Oe/s and 200 Oe/s 
are shown in Fig.ll. At 2 K, the SP effect is obviously ob- 
served, and by increasing temperature, the second-peak 
shifts to the low magnetic field. The very sharp mag- 
netization peak was also found near the zero field. The 
symmetric curves suggest that the bulk pinning instead 
of the surface barrier dominates in the sample. 

Based on the Bean critical state modelfSS, it is known 
that the superconducting transient current js oc AM. 
The non-monotonic field and temperature dependence 
of is was clearly observed in the SP region (shown in 
Fig. 12). At a low field, magnetization relaxation rate 
has an upturn when decreasing the magnetic field, which 
is in accord with the valley of the MHL. In the inter- 
mediate field region, Q decreases with the increasing of 



T 




H(T) 

FIG. 12: (Color online) Magnetic field dependence of the tran- 
sient superconductivity current density is defer from the width 
AAl and dynamic magnetization relaxation rate from 2 K to 
16 K with the sample x= 0.12. 

the magnetic field, this indicates a crossover to a vor- 
tex glass region with relatively strong pinning of vortex. 
Above the second-peak maximum value Hsp, the value of 
Q increases drastically showing a crossover to the regime 
of plastic motion. 

The temperature dependence of Q and js for the over- 
doped sample x = 0.12 was presented in Fig. 13. Obvi- 
ously, there is still a plateau in the intermediate temper- 
ature region for each Q — T curve, which is corresponding 
to the linear region of log — T curve. And the region of 
the plateau shrinks quickly with increasing the magnetic 
field. This behavior can be understood within the picture 
of collective vortex pinning model. In the high tempera- 
ture region, the Q increases abruptly, and the log js drops 
down drastically, which again shows a crossover from the 
collective elastic motion to plastic motion. 

The similar vortex behavior with the optimally doped 
sample can be derived from the vortex phase diagram for 
the sample with x = 0.12 (shown in Fig. 14). From 
the magnetic measurements, three characteristic fields 
are determined as shown by the filled symbols in Fig. 
14. The SP field Hsp locating at the peak of the js — H 
curve shown as the solid pentagram now is much lower 
than that in the optimally doped sample. The irre- 
versibility field Hirr is determined by taking a criterion 
of 0.1 emu/cm'^ shown as the solid star. The upper crit- 
ical field Hc2 with 95 percent p„ is shown as the filled 
circles. The Hsp — T, Hirr — T curves are well fitted 
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FIG. 13: (Color online) Temperature dependence of log and 
the dynamic magnetization relaxation rate Q with x — 0.12 
at 0.5, 1, 2, 4 and 6 T. 
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T(K) 



FIG. 14: (Color online) The vortex phase diagram of the 
overdoped sample {x — 0.12). Three characteristic fields Hsp, 
Hirr and Hc2 are shown. 



by the expressions Hsp{T) = Hsp[0){l - T jT^f-^ and 
Hirr{T) = i?irr(0)(l - T/Tcf '^. Below the Hsp Hne the 
vortex behavior is consistent with the collective creep 
model. Between H^p line and Hi„ line, it can be under- 
stood by the model of plastic motion. Above the i/j^r 
line, the vortex liquid phase exists. 
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H(T) 



FIG. 15: (Color online) The MHLs of underdoped sample (a; 
= 0.06) from 2 K to 10 K. The solid line represents the data 
measured at the magnetic sweeping rate of 200 Oe/s while 
the dotted hne at 50 Oe/s. 



C. Underdoped sample with x = 0.06 

Fig. 15 shows the magnetization hysteresis loops of the 
sample x — 0.06 measured at 2-10 K with the magnetic 
field sweeping rate 50 Oe/s and 200 Oe/s, respectively. 
At all temperatures we measured here, the SP effect can- 
not be obviously observed. As the optimal doped and 
overdoped sample the very sharp magnetization peak 
was also found near the zero field. The symmetric MHL 
curves suggest that the bulk pinning instead of the sur- 
face barrier dominates even in this very underdoped sam- 
ple. 

The field dependence of js and Q from 2 K to 10 K was 
presented in Fig. 16. The js decreases monotonously with 
the increase of magnetic field, without showing a SP ef- 
fect. Correspondingly, Q increases with the increasing of 
the magnetic field rapidly and quickly reaches the upper 
limit 100%. This indicates that the vortex motion in the 
most temperature regime investigated for this sample is 
dominated by plastic motion. The temperature depen- 
dence of js and Q measured at different magnetic fields 
were presented in Fig. 17. By increasing temperature, ini- 
tially log js decrease linearly, while it quickly evolves into 
a quick drooping down in the high temperature region. 
In the low field region (0.5 T and 1 T), the relaxation rate 
Q rises relatively slowly with temperature in low and in- 
termediate temperature region, while at 2 T, 3 T, 4 T 
and 5 T, the Q curve increases dramatically on warm- 
ing, which suggests the quick evolvcmcnt of the plastic 
motion of vortices. 

The vortex phase diagram of the underdoped sample 
with X — 0.06 is presented in Fig. 18. Between the upper 
critical field and the irreversibility line, it is the region 
for vortex liquid, while there is now a region below the 
irreversible line which may have a low dissipation, but be 
dominated by plastic motion. It is still unclear whether 
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H(T) 

FIG. 16: (Color online) Magnetic field dependence of the 
MHL width AM and dynamic magnetization relaxation rate 
measured for the underdoped sample a;=0.06 from 2 K to 10 
K. 



in this sample we can still see the SP effect at tempera- 
tures below 2 K. It is obvious that the SP effect is absent 
at temperatures above 2 K. One of the interpretations is 
that the elastic energy which is based on the characteris- 
tic pinning energy is too weak to sustain an elastic object 
in the vortex system. Therefore in this sample, the vortex 
motion below Hirr(T) is dominated by a plastic manner. 
While we should not exclude the possibility that in the 
regime at much lower temperatures, the elastic motion 
still exists. 



IV. DISCUSSION 

In all six samples investigated here, sizable magnetiza- 
tion relaxation rate has been observed, which indicates 
a relatively small characteristic pinning energy. In the 
samples around the optimal doping level {x « 0.08) a 
SP effect was easily observed, while it becomes invisible 
in the very underdoped samples, and hardly visible for 
the highly overdoped one. The missing of the SP effect 
in the very underdoped and overdoped regions is under- 
standable since the pining energy becomes much weaker 
than that of the optimally doped sample. In this case, the 
dislocation in the vortex system can easily move leading 
to a plastic motion as appeared in the high field region of 
the optimally doped sample. Therefore we believe that 
the missing of the SP effect in the very underdoped or 
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FIG. 17: (Color online) Temperature dependence of logjs 
and dynamic magnetization relaxation rate Q with x — 0.06 
at 0.5, 1, 2, 3, 4 and 5 T. It is clear that the plateau of the 
relaxation rate becomes very short followed by a quick rising 
up of relaxation rate. 
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FIG. 18: (Color online) The phase diagram of the ounder- 
doped sample (a;=0.06) 

overdoped samples is induced by the missing of the elas- 
tic collective pinning and creep. The monotonic rising of 
the relaxation rate in whole temperature region for the 
very underdoped and overdoped samples support this ar- 
gument. Another interesting observation is that a very 
sharp magnetization peak was observed near the zero 
field in all samples which is found to correspond very well 
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FIG. 19: (Color online) Temperature dependence of (a) the 
log j's, Q and T /Q for all six samples investigated in this work. 
In low temperature region, a roughly linear temperature de- 
pendence of logj3(T) vs. T can be observed for all samples 
besides the one with x = 0.06. 



to a much reduced relaxation rate. Fig. 19 shows the tem- 
perature dependence of Q, js and T/Q with all doping 
levels at 0.5 T. We found an bell-shaped or a plateau-like 
temperature dependence of Q with moderate relaxation 
rate in most samples, which can be understood within 
the collective creep theory. The plateau region in the 
very underdoped x = 0.06 and very overdoped sample 
X — 0.15 is very short because of the low pinning en- 
ergy. According to eq.6, the value of T/Q in the T — 
K approach would give the characteristic pinning energy 
C/c(0). Fig. 20 shows the doping dependence of ?7c(0) at 
0.5 T with all doping levels. The linear logjs-T curves 
of all doped samples at 0.5 T suggest the thermally acti- 
vated collective vortex creep feature in the intermediate 
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FIG. 20: (Color online) Doping dependence of jco, Uco and 
Tc for all six samples investigated here. It is clear that both 
jco and UcO have a very similar doping dependence of Tc. 
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FIG. 21: (Color online) Doping dependence of the Hirr and 
Hc2 for all six samples. 
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temperature region. The jc(0), Uc{0) and Tc curves have 
the similar dome shape with doping, this may suggest 
that the vortex pinning is through the disorders with 
small condensation energy. By taking a GIS treatment 
on the magnetization data, we found that the vortex pin- 
ning is probably achieved through the spatial fluctuation 
of the superconducting transition temperatures, the so- 
called (JTc— pinning. Further investigations are strongly 
recommended to clarify this issue. With the decreasing of 
Tc, the samples have lower irreversible field and the upper 
critical field, as shown in Fig 21. While compared to the 
underdoped samples, the overdoped samples with similar 
Tc seem to have higher irreversible field and the upper 
critical field. The SP in underdoped samples disappear 
more quickly with doping than that in overdoped sample. 
For example, when the doping level drops down from x 
— 0.07 to 0.06, the SP effect disappears, while it is still 
visible when the doping level goes to 0.12, although very 
weak. This tendency is actually similar to the doping 
dependence of Tc and the characteristic pinning energy 
Uc{0), as shown in Fig. 20. This again suggests that 
the SP effect is dependent on the pinning energy which 
governs actually the threshold of the plastic motion of 
vortex system. 

V. CONCLUDING REMARKS 

We measured magnetization and its relaxation in 
Ba(Fei_a;Coa;)2As2 single crystals at various doping lev- 
els ranging from very underdoped to very overdoped 
regime. Detailed analysis lead to the following major 
conclusions: 

(1) In all samples, sizable magnetization relaxation 
rate has been observed, which suggests relatively weak 
vortex pinning energy. The characteristic pinning energy 
obtained here for the optimally doped sample is about 



100 K at about 0.5 T. 

(2) A very sharp magnetization peak was observed near 
zero field which corresponds to a much reduced relaxation 
rate. This may be induced by the extremely non-linear 
B{x) near the edge when the external field is swept to 
zero, or due to the surface barrier for the vortex exit and 
entering at the edge. 

(3) The second-peak effect was easily observed in the 
samples around the optimal doping level {x «0.08), but 
it becomes hardly visible in the very underdoped and 
highly overdoped samples. We attribute the missing of 
the SP effect to the much weaker pinning energy, which 
leads to a plastic motion of vortices in wide temperature 
regions. Through the SP region the transient supercon- 
ducting current density shows the non-monotonic field 
and temperature dependence. 

(4) The weak temperature dependence of relaxation 
rate together with the treatment of the Generalized- 
Inversion-Scheme, point to the fact that the model of 
collective vortex pinning and creep works well in describ- 
ing the vortex dynamics in iron pnictides. The vortex 
pinning is probably achieved through the spatial fluctu- 
ation of the transition temperatures, which would mean 
an intrinsic inhomogeneity of the iron-pnictide supercon- 
ductors. 
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